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Abstract. The theory and phenomenology of light sterile neutrinos at the eV mass 
scale is reviewed. The reactor, Gallium and LSND anomalies are briefly described and 
interpreted as indications of the existence of short-baseline oscillations which require 
the existence of light sterile neutrinos. The global fits of short-baseline oscillation data 
in 3+1 and 3+2 schemes are discussed, together with the implications for /3-decay 
and neutrinoless double-/3 decay. The cosmological effects of light sterile neutrinos are 
briefly reviewed and the implications of existing cosmological data are discussed. The 
review concludes with a summary of future perspectives. 
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1. Introduction 

The possible existence of sterile neutrinos [l], which do not have the standard weak 
interactions of the three active neutrinos z+, and u T1 is currently a hot topic of 
theoretical and experimental research which could provide valuable information on the 
physics beyond the Standard Model (see Refs. 12 20]). 

Sterile neutrinos are singlets of the Standard Model gauge symmetries which can 
couple to the active neutrinos through the Lagrangian mass term. In practice there 
are bounds on the active-sterile mixing, but there is no bound on the number of 
sterile neutrinos and on their mass scales. Therefore the existence of sterile neutrinos 
is investigated at different mass scales. This review is devoted to the discussion of 

1 This review is dedicated to the memory of Hai-Wei Long, our dear friend and collaborator, who 
passed away on 29 May 2015. He was an exceptionally kind person and an enthusiastic physicist. We 
deeply miss him. 
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sterile neutrinos at the eV scale, which can explain the anomalies found in some short- 
baseline neutrino oscillation experiments (discussed in details in Section [3j) . However, 
it is important to remind that there are other very interesting possibilities which are 
under study: very light sterile neutrinos at a mass scale smaller than 0.1 eV, which could 


affect the oscillations of solar 21-23 and reactor 24-28 neutrinos; sterile neutrinos at 


the keV scale, which could constitute warm dark matter according to the Neutrino 
Minimal Standard Model (z/MSM) |29]-[33| (see also the reviews in Refs. [8, 9,|T8,[T9]); 


sterile neutrinos at the electroweak scale |34,35 or above it 35,36 , whose effects may 


be seen at LHC and other high-energy colliders. Let us also note that there are several 


interesting models with sterile neutrinos at different mass scales 37--52 

The plan of this review is as follows. In Section [2] we introduce the main 
theoretical aspects of neutrino mixing and oscillations in the extensions of standard 
three-neutrino (3u) mixing with light sterile neutrinos. In Section [3] we summarize 
the three experimental indications in favor of short-baseline neutrino oscillations which 
require the existence of sterile neutrinos: the reactor, Gallium and LSND anomalies. 
In Section [4] we discuss the global fits of short-baseline neutrino oscillation data in the 
framework of the simplest mixing schemes with one or two sterile neutrinos at the eV 
scale. In Sections [5] and [6] we review the effects of eV-scale massive neutrinos in /? decay 
and neutrinoless double-/? decay and in cosmology. Finally, in Section [7] we present our 
conclusions and our view of the perspectives of the search for eV-scale sterile neutrinos. 


2. Neutrino mixing and short-baseline oscillations 


The Standard Model of electroweak interactions 53 -55 based on the SU(2 )l x U(l)y 
gauge symmetry is a superb theory which can explain the majority of terrestrial 
experimental observations. However, it does not account for neutrino masses, whose 
existence have been proved without doubt by the measurement of neutrino oscillations 
in solar, atmospheric and long-baseline neutrino oscillation experiments (see Refs. 17, 
56 60]). The simplest way to extend the Standard Model in order to take into account 


neutrino masses is through the introduction of SU(2) L x U(l)y singlet fields which 
are traditionally called “right-handed neutrino” fields or “sterile neutrino” fields. The 
adjective “right-handed” indicates that they do not belong to SU(2)i left-handed 
multiplets. Therefore, they are “sterile”, because they do not have Standard Model 
weak interactions. Moreover, assuming that they have zero hypercharg^J they are 
neutral and can be called “neutrino” fields. Many models which extend the Standard 
Model include these right-handed sterile neutrino fields (see Refs. (3[ 8,12,18; 63, 64 ). 


In the following we consider the general theory of neutrino mixing in which we have 
the three standard active left-handed flavor neutrino fields u eR , u t l and N s sterile 
right-handed flavor neutrino fields is slR , ..., n NaR . The most general Lagrangian mass 
term which can be written with these fields is (the superscript “(F)” indicates the flavor 

2 We do not consider here the exotic possibility of a small nonzero hypercharge of the right-handed 
neutrino fields, which would imply that neutrinos are Dirac and millicharged particles (see Refs. 61 62 
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C nmss ^ 2 i/p T CKMu ( P + R.c., ( 1 ) 

where (the superscripts “(a)” and “(s)” indicate, respectively, the column matrices of 
active and sterile neutrino fields) 



( a ) 

M = 


(Vel} 
VrL) 


(s) c 

V R = 


'si R 


\ 


‘ U s Na R) 


( 2 ) 


and C is the unitary charge-conjugation matrixrl such that C 7 J C 1 = — 7 ^ and 

' - - rj 1 

C 1 = —C. For any held ip the charge-conjugated held -0 C is given by -0 C = Cip and 
charge conjugation transforms the chirality of a held (e.g. i\) c R is left-handed). In general, 
M is a complex symmetric mass matrix, which can be diagonalized with the unitary 
transformation (the superscript “(M)” indicates the mass basis) 




with 




( 3 ) 


where N = 3 + N s is the total number of neutrino helds. The matrix U is a N x N 
unitary matrix such that 


U 1 MU = diag(mi,..., mjv), 


( 4 ) 


with real and positive masses (see Refs. [56,65]). The Lagrangian mass 

term (|TJ) becomes 


i N 1 N _ 

jO-mass = - m k^lrP 1J kL + H.c. = - - ^ m k u^ L v kL + H. 


N 




- E 


rn k v k v k , (5) 


k =1 


k =1 


k =1 


with the massive Majorana neutrino helds v k = v kkj + u < j tL which satisfy the Majorana 
constraint v k = Hence, in the general case of active-sterile neutrino mixing the 
massive neutrinos are Majorana particletQ 

The physical effects of the unitary transformation ([3]) are due to the non-invariance 
of the weak interaction Lagrangian. Let us hrst consider the leptonic charged-current 
weak interaction Lagrangian. In the havor basis where the mass matrix of the charged 
leptons i e = e, = /i, l T = r, is diagonal, we have 


N 


Ccc = ~7 5 X) ZolYvolWI + H.c. = t aL l p U ak v kL Wl + H.C.. 

v a=e,ii,T v fc=l a=e,ii,T 


( 6 ) 


56 


3 We use the notations and conventions in Ref. 

4 However, it is not excluded that the mixing is such that there are pairs of Majorana neutrino fields 
with exactly the same mass which form Dirac neutrino fields. 
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£°° = -^=4-Ai‘ ) W'2 + h.c. = + H.c., (7) 


with 



4 a) = I/4 M) and U = U\ 3xN . 


( 8 ) 


V/ 

The mixing matrix U is a 3 x At rectangular matrix formed by the first three rows of 
hi. Therefore, the number of physical mixing parameters is smaller than the number 
necessary to parameterize the unitary matrix hi. This is due to the arbitrariness of the 
mixing in the sterile sector, which does not affect weak interactions. A careful analysis 
(see Ref. [56j) shows that the 3 x N mixing matrix U can be parameterized in terms of 
3 + 3N S mixing angles and 3 + 3 N s physical phases, of which 1 + 2N S are Dirac phases 
and N — 1 are Majorana phases. For such parameterization, it is convenient to use the 
scheme 


U 


N 


II ll Wab / ' i3 " l3/ ' 12 


ia=l b=4 


diag(l, 


21 


,e 


i A a 


J 3 xN 


( 9 ) 


The unitary N x N matrix W ab = W ab (9 abl ijab) represents a complex rotation in the a-b 
plane by a mixing angle 9 ab and a Dirac phase rj ab . Its components are 


[W ab (‘9 ab , Vab)] rs = $rs + (c ab - 1) ( 5 r J sa + 8 rb S sb ) +s ab (e lVab 8 ra 5 sb - e tVab 8 rb 8 sa ) , (10) 


where c ab = cos i9 ab and s a b = sin 'd ab . The order of the product of W ab matrices in 
Eq. (J 9 ]) is arbitrary. The orthogonal matrix R ab = W ab (9 ab ,0) represents a real rotation 
in the a-b plane. The square brackets with subscript 3 x N indicate that the enclosed 
N x N matrix is truncated to the first three rows. The Majorana phases A 21 ,... Ajvi., 
which are physical only if massive neutrinos are Majorana particles, are collected in a 
diagonal matrix on the right0 Moreover, not all the phases r] ab in the product of W ab 
matrices in Eq. (J9]) are physical, but one can eliminate an unphysical phase for each 
value of the index 6 = 4,..., iV (see Ref. [561). 

The scheme ([9]) has the advantage that in the limit of vanishing active-sterile mixing 
the mixing matrix reduces to the three-neutrino (3zz) mixing matrix in the standard 
parameterization 


U [ 3u) = [R 23 W 13 R 12 ] 3x3 diag(l,e iA21 ,e a31 ) 

^ C12D3 S12C13 

-S12C23 - Ci2S23Sl3e l,?13 C12C23 - Si 2 S23Sl 3 e lrl13 

\ 5 i 2 s 23 - C 12 c 23 s 13 e ir113 -C12S23 - Si 2 c 23 Si3e im3 


s 13 e~ im3 \ 

S23C13 

C23 C 13 ) 


(l 

0 

v° 


0 

0*^21 


0 


0 ^ 
0 


e* Asi / 

( 11 ) 


5 It is possible to choose any other diagonal matrix with N — 1 phases, as for example 
diag(e lAl2 , l,e lA32 ,... ,e* AjV2 ), etc. 
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It is convenient to choose in Eq. ([9]) the order of the real or complex rotations for 
each index 6 > 4 such that the rotations in the 3-6, 2-6 and 1-6 planes are ordered 
from left to right. In this way, the first two lines, which are relevant for the study of 
the oscillations of the experimentally more accessible flavor neutrinos v e and are 
independent of the mixing angles and Dirac phases corresponding to the rotations in all 
the 3-6 planes for 6 > 4. Moreover, the first line, which is relevant for the study of v e 
disappearance, is independent also of the mixing angles and Dirac phases corresponding 
to the rotations in the 2-6 planes for 6 > 3. For example, one can choose 

U = [W 3N R 2N W 1N ■ ■ ■ W U R 24 W U R 23 W 13 R 12 ] 3xN diag(l, e iAzi ,... ,e iAjvl ) , (12) 


or 


U = [W' 


3N ... w u W 2N ■ ■ ■ W 24 R 1N 


R u R 23 W 13 R 12 ] 3xN diag(l, e* Azi ,... ,e Ujvl ) . 

(13) 

Let us now consider the neutrino neutral-current Lagrangian 


£nc — 
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(a) p (a) ry 
-VlY^l Z P = 




(14) 


2cos'$w L L p 2cos'$w L 

Since the rectangular 3 x N mixing matrix U is formed by the first three rows of the 
unitary matrix U , we have 


Utf = 1 3 x 3 , but YU ^ l NxN . 


(15) 


Therefore, the GIM mechanism |66| does not work in neutral-current weak interactions 
[671 and it is possible to have neutral-current transitions among different massive 
neutrino^ 

The introduction of sterile neutrinos is allowed by the fact that it has no effect or 
small effect^] on the effective number of active neutrinos which contributes to the decay 
of the Z-boson. This number has been determined with high precision to be close to 
three by the LEP experiments |72|: 


NY = 2.9840 ± 0.0082. 


(16) 


(Z) 

In this review we will consider sterile neutrinos at the eV scale, for which N„ ; is given 

(17) 


by 70,71 


N 


N m = £ 

j,k=l 


= 3. 


E u * u <* 

a=e,/i,T 

Hence, there is no constraint on the number and mixing of these light sterile neutrinos 
from the high-precision LEP measurement of N„ . 

6 This is a special case of the general theorem that the weak leptonic neutral current is nondiagonal 
in the mass basis if the leptons of a given charge and chirality have different weak isospins 68 . 

7 Sterile neutrinos at mass scales larger than the muon mass affect the determination of the Fermi 


constant Gp through muon decay 69 . Those at mass scales larger than 771^/2 can induce a kinematical 
suppression of nY 70,71 . 













Light sterile neutrinos 


6 


parameter 

mass 

order 

best 

fit 

la : 

range 

2a : 

range 

3a 

range 

Ami OL /10- 5 eV 2 


7.54 

7.32 

- 7.80 

7.15 

- 8.00 

6.99 

- 8.18 

sin 2 t?i 2 /10” 1 


3.08 

2.91 

- 3.25 

2.75 

- 3.42 

2.59 

- 3.59 


NO 

2.43 

2.37 

- 2.49 

2.30 

- 2.55 

2.23 

- 2.61 

Am ATM/ 10 ~ 3 ey2 

IO 

2.38 

2.32 

- 2.44 

2.25 

- 2.50 

2.19 

- 2.56 


NO 

4.37 

4.14 

-4.70 

3.93 

-5.52 

3.74 

- 6.26 

sin 2 '^ 23/10 1 

IO 

4.55 

4.24 

- 5.94 

4.00 

- 6.20 

3.80 

- 6.41 


NO 

2.34 

2.15 

- 2.54 

1.95 

- 2.74 

1.76 

- 2.95 

sin 2 r? 13 /10 2 

IO 

2.40 

2.18 

- 2.59 

1.98 

- 2.79 

1.78 

- 2.98 


Table 1 . Values of the neutrino mixing parameters obtained in Ref. 73 with a 
global analysis of neutrino oscillation data in the framework of three-neutrino mixing 
with the normal ordering (NO) and the inverted ordering (10). 


The measurements of neutrino oscillations in solar, reactor and accelerator 
experiments determined the existence of two squared-mass differences: the solar and 
atmospheric squared-mass differences 

AmgoL — 7.5 x lCT 5 eV 2 , Am 2 ATM ~ 2.4 x 10 _ 3 eV 2 . (18) 

It is convenient to label the three light neutrino masses according to the convention 

a "4 ol = A m 2 21 < A m 2 ATM = i | A m 2 31 + Am 2 32 \ , (19) 

with A m 2 - k = nij — m 2 k . The absolute value in the definition of A m\ TM is necessary, 
because there are two possible orderings of the neutrino masses: the normal ordering 
(NO) with mi < m 2 < m 3 and Am| 1? A m 32 > 0; the inverted ordering (10) with 
m 3 < mi < m 2 and A m 3ll A m \ 2 < 0 . 

Table [ljshows the results of the determination of the mixing parameters obtained 
in Ref. |73| from a global fit of neutrino oscillation data (see also Refs. |74,75]). The 
largest uncertainty is that of -$ 23 , which is known to be close to maximal (7t/4), but it is 
not known if it is smaller or larger than 7t/4. For the Dirac CP-violating phase 7733 , there 
is an indication in favor of 7/13 ~ 37t/2, which would give maximal CP violation, but at 
3 cr all the values of 7713 are allowed, including the CP-conserving values 7713 = 0, 7 r. 

The standard framework of mixing can be extended with the introduction of non¬ 
standard massive neutrinos only if their mixing with the active neutrinos is sufficiently 
small in order not to spoil the successful 3z/ mixing explanation of solar, atmospheric 
and long-baseline neutrino oscillation measurements discussed above. In other words, 
the non-standard massive neutrinos must be mostly sterile, i.e. 

\U ak \ 2 « 1 (a — e, fi,T] k — A,..., N). (20) 

In the following we will always assume this constraint. 
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In this review we consider mainly the so-called 3+1 scheme in which there is a 
non-standard massive neutrino (mostly sterile) at the eV scale which generates a new 
squared-mass difference 


Am 


SBL 


1 eV , 


( 21 ) 


in order to explain the anomalies found in some short-baseline (SBL) neutrino oscillation 
experiments (see Section [3]). We assume that the three standard massive neutrinos are 
much lighter than the eV scale. We will consider also the so-called 3+2 scheme in 
which there are two non-standard massive neutrinos (mostly sterile) at the eV scale 
and the 3+1+1 scheme in which there is a non-standard massive neutrino at the eV 
scale and another at a larger scale. We do not consider schemes in which Amg BL is 
obtained with one or more very light (or massless) non-standard massive neutrinos and 
the three standard massive neutrinos have almost degenerate masses at the eV scale (e.g., 
the 1+3, 1+3+1 and 2+3 schemes), because this possibility is strongly disfavored by 
cosmological measurements 1 76 1 and by the experimental bound on neutrinoless double-/? 
decay (assuming that massive neutrinos are Majorana particles; see Ref. |77|). Figure [l] 
shows a schematic illustration of the 3+1, 3+2 and 3+1+1 neutrino mixing schemes 
taking into account for each scheme the two possible mass orderings of the three lightest 
standard neutrinos. Let us emphasize that these mixing schemes must be considered 
as effective, in the sense that the existence of more non-standard massive neutrinos is 
allowed, as long as their mixing with the three active neutrinos is sufficiently small to 
be negligible in the analysis of the data of current experiments. 

For the study of neutrino oscillations in vacuum it is convenient to use the following 
general expression of the probability of l z/„ —> < vg oscillations 


78 


79): 


R 


<-) (-) 


<W - 4^ \U ak \ 2 (5 a p - \Up k \ 2 ) sin 2 A 


kp 


k^p 


(+) 


+ 8 ^ ^ | UajUfijUakUpk | sin Afcp sin Ajp cos( Aj k rjnrijk) • 


( 22 ) 


j>k 

ik+v 


where 


^kp 


Am\pL 


4 E 


VafSjk = arg [K j U f > j U ak Up k ] , 


(23) 


and p is an arbitrary fixed index, which can be chosen in the most convenient way 
depending on the case under consideration. In the case of three-neutrino mixing, there 


is only one interference term in Eq. (22), because for any choice of p there is only one 


possibility for j and k such that j > k. 

We are interested in the effective oscillation probabilities in short-baseline 
experiments, for which A 21 -C A 31 -C 1. Let us consider the general 3 +N S case in 
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-- Vl 
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- Vl 


A m A TM 
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Am | 4 
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Figure 1. Schematic illustration of the 3+1, 3+2 and 3+1+1 neutrino mixing 
schemes taking into account for each scheme the two possible mass ordering, normal 
(NO) and inverted (10), of the three lightest standard neutrinos. In the two 3+2 
schemes Amiq ss Amg BL , whereas in the two 3+1+1 schemes Am^ Amg BL . 


which ~ Am| BL and A/, :1 rs 1 for k > 4. Choosing p — 1 in Eq. (23), we obtain 


N 


p(SBL) _ 

(-) (-) — 
Va^-Vp 


S a p - \U»k\ 2 (<W - \Upk\ 2 ) sin 2 Afci 


fc=4 
N N 


(+) 


\UajUpjU ak Up k \ sin Afci sin A^i cos(A jfc - 7 ) a/3jk ). (24) 

k=4 j=k-{-1 

Considering the survival probabilities of active neutrinos, let us define the effective 
amplitudes 


sin" 


= 4|C/„ t | 2 (1 - \U ak \ 2 ) ~ i\U ak \ 2 (a = e,p.,r, k > 4), 


(25) 


where we have taken into account the constraint in Eq. (20). Dropping the quadratically 


suppressed terms also in the survival probabilities, we obtain 

N 

^ „:„2 o -. offc ') „:„2 a 

'kl 


f 2) - 1 - X] sin2 2 ^a sin 2 A a 

Vol ^Vcti k =4 


(a = e, p, t). 


(26) 


Hence, each effective mixing angle $aa parameterizes the disappearance of'u'a due to its 
mixing with v k . 

Let us now consider the probabilities of short-baseline l v a transitions between 

two different active neutrinos or an active and a sterile neutrino. We define the transition 
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' ah 


! | u pk 


amplitudes 

sin 2 20$ = 4| U t 
which allow us to write the transition probabilities as 

N 


M /3; k> 4), 


(27) 


p(SBL) _ 

(-) (-) — 
Vc^-lSp 


^ sin 2 20$ sin 2 A kl 


k =4 


iC?) 


(+) 


AT TV 

2 EE sin 20$ sin 20$ sin A fci sin A j4 cos(A jfc - r] a/3jk ). (28) 

/c=4 jf=fc-f-l 


From the first line one can see that each effective mixing angle 0$ parameterizes the 
amount of*//^ —» ( i/g transitions due to the mixing of^ and ( zyg witl/iy,. The second line 
in Eq. (28) is the interference between the l u k and V] contributions, which depends on 
the same effective mixing angles. 

Considering now the transitions between two different active neutrinos, from 
Eqs. (25) and (27) one can see that for each value of k > 4 the transition amplitude 
sin 20$^ and the disappearance amplitudes sin 2 -dal and sin 20$ depend only on the 
elements in k th column of the mixing matrix and are related bj 


82 


sin 2 20$ ~ - sin 2 20gJ sin 2 20$ 


(a = e,p,r). 


(29) 


This relation is very important, because it constrains the oscillation signals that can 
be observed in short-baseline appearance and disappearance experiments in any 3 +N S 
mixing scheme with sterile neutrinos. Its experimental test is crucial for the acceptance 
or rejection of these schemes. In particular, since both sin" 20+i and sin 20$ are small 
for a, (3 = e, /j,,t the amplitudes of the short-baseline transition probabilities between 
active neutrinos are quadratically suppressed. We will see in Section [4] that the current 
short-baseline data have an appearance-disappearance tension due to the constraint in 
Eq. 

In the following part of this Section we discuss briefly the main peculiar 
characteristics of short-baseline oscillations in the cases of 3+1, 3+2 and 3+1+1 neutrino 
mixing. 

3+1 mixing. In the case of 3+1 neutrino mixing [5j80,8T|j83|, we have Am^ = Am| BL 
and A 41 ~ 1 in short-baseline experiments. The transition and survival probabilities 
can be written as 


^ sin 2 2 d a 0 sin 2 A 4i (a^), 

Va-yt/f 3 


~ 1 - sin 2 2 i) c 


, sin 2 A 


41, 


with the transition and survival amplitudes 

sin 2 — 41C/^41 2 11 2 sin 2 2tf„ = 4|t/ t>4 | 2 (1 - \U aa \ 2 ) , 


(30) 


( 31 ) 


This relation was derived in the case of 3+1 mixing (see Eq. (32)) in Refs. 80 


81 
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and with the appearance-disappearance constraint 80,81 


sin 2 2 d a j 3 ~ - sin 2 2 d aa sin 2 2 , dpp (a = e, p, r). 


(32) 


The transition and survival probabilities in Eq. (30) depend only on the largest 


squared-mass difference = Amg BL and on the absolute values of the elements in 

the fourth column of the mixing matrix. The transition probabilities of neutrinos and 
antineutrinos are equal, because the absolute values of the elements in the fourth column 
of the mixing matrix do not depend on the CP-violating phases in the mixing matrix. 
Hence, even if there are CP-violating phases in the mixing matrix, CP violation cannot 
be measured in short-baseline experiments. In order to measure the effects of these 
phases it is necessary to perform experiments sensitive to the oscillations generated by 
the smaller squared-mass differences Am} TM 


84-86 


or A m 2 SOh 


87 


3+2 mixing. In the case of 3+2 neutrino mixing 88-92 , we have Am? 


"51 


Am 2 ! = 


AmgBL an d A r 


*51 ~ A41 ~ 1 in short-baseline experiments. From Eq. (26) we obtain the 
short-baseline survival probabilities of active neutrinos 

— 1 - sin 2 2$£ l 2 sin 2 A 41 - sin 2 2d£2 sin 2 A 51 (a = e, p, r), (33) 

Va^fVa 


and from Eq. (28) we obtain the short-baseline transition probabilities 


P ( { ! f J2) - sin 2 2d^ sin 2 A 4i + sin 2 2d 


j(4) 


sin 2 A 


51 




+ 2 sin 2^1 sin 2$^ sin A 41 sin A 51 cos(A 54 - 770,54) (a ^ j8). (34) 

The appearance and disappearance amplitudes are related by the general constraint 


in Eq. (29). The 3+2 scheme has the important characteristic that CP violation is 


observable in short-baseline experiments through the asymmetries 

^(sm) = p ( smq _ p (s Biq _ 4 gin 2 ,(4) gin 2 ,(g gin Aji sin a 51 sin A 54 sin r] af}54 , (35) 

for a 7^ /3. 


3+1+1 mixing. I 11 the case of 3+1+1 mixing 93-96 , we have Am 2 x Am| x = 
51 )+ A 4 i ~ 1 in short-baseline experiments. The corresponding oscillation 


AmgBL an d A 


probabilities can be obtained from those in the case of 3+2 mixing by averaging the 


oscillations due to Am^: 


++ — 1 — sin 2 sin 2 A 4 , - - sin 2 2i(a = e, /!, r). 


(36) 


and 


P 


(SBL) 


(-) (-) 
V P 


~ sin 2 2d^j sin 2 A 4i -|— sin 2 2d 


(5) 

a/3 


+ sin 2 ^ sin 2 d^ sin A 4i sin(A 4i V 77 ^ 54 ) (a ^ P). (37) 

Hence, in the analysis of short-baseline data in the 3+1+1 scheme there is one effective 
parameter less than in the 3+2 scheme (Am|i), but CP violation effects generated by 
the phases 7p, 54 are observable. 
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Figure 2. Ratios R of the measured (iV exp ) and calculated (JV ca i) number of electron 
antineutrino events in reactor experiments at different distances L. The horizontal 
shadowed red band shows the average ratio R and its uncertainty. For each experiment 
the error bar shows the experimental uncertainty. The values of the ratios of the 
long-baseline experiments Daya Bay, Double Chooz, Chooz and Palo Verde have been 
obtained by subtracting the effect of $ 13 -driven oscillations. 


3. Short-baseline anomalies and constraints 


In this Section we review the three experimental indications in favor of short-baseline 
neutrino oscillations, which require the existence of at least one additional squared- 
mass difference, Amg BL , which is much larger than Arog OL and A rn\ TM \ the reactor 


antineutrino anomaly in Subsection 3.1, the Gallium neutrino anomaly in Subsection 3.2 


and the LSND anomaly in Subsection 3.3 


3.1. The reactor antineutrino anomaly 


The reactor antineutrino anomaly 97 , is a deficit of the rate of u e observed in several 


short-baseline reactor neutrino experiments in comparison with that expected from the 
calculation of the reactor neutrino fluxes |12,|98 ; ! , 99]. It has been analyzed by several 


authors 12,27,97,100 111 


The statistical significance of the anomaly depends on the estimated uncertainties 
of the calculated reactor antineutrino fluxes, on which there have been some debate 


112-119 , especially after the discovery of an excess at about 5 MeV of the reactor 


antineutrino spectrum in the RENO 120 , Double Chooz 121 and Daya Bay 122-124 


experiments. However, since this problem is controversial and far from being solved, in 
the following we assume the values and uncertainties of the reactor antineutrino fluxes 
presented in Refs. 12,98,99 , which have been obtained from the available database 


information on nuclear decays and from the electron spectra associated with the fission 
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51 Cr 

37 Ar 

E [keV] 

747 

752 427 

432 

811 813 

branching ratio 

0.8163 

0.0849 0.0895 

0.0093 

0.902 0.098 


Table 2. Energy (E) and branching ratio of the neutrino lines produced in the 
electron-capture decays of 51 Cr and 37 Ar. 


of 235 U, 239 Pu, and 241 Pu measured at ILL in the 80’s 12511128 . The electron spectra 


239i 


associated with the fission of 235 U and 238 U has been measured recently at the scientific 
neutron source FRM II in Garching [l291. In this experiment the 235 U electron spectrum 


has been normalized to that measured at ILL 126,128 in order to reduce the systematic 
uncertainties. The 238 U antineutrino spectrum obtained from the conversion has only 
a spectral distortion of about 10% with respect to the standard 238 U antineutrino 
spectrum (12.08.!)!)]. 

In reactor neutrino experiments electron antineutrinos are detected through the 
inverse neutron decay process 

v e + P ~> n + e + (38) 

in liquid-scintillator detectors. The cross section of the detection process is ap eP (E e ) oc 
E e p e (see Refs. 56,130, 13lj ), where E e and p e are, respectively, the positron energy 
and momentum. Neglecting the small recoil energy of the neutron, the neutrino energy 
E can be calculated from the measurable kinetic energy T e of the positron through the 
relation 

E ~ T e + m e + m n — rn p ~ T e + 1.8 MeV, (39) 

where rn p and m n are, respectively, the proton and neutron masses. Hence, the neutrino 
energy threshold for the detection process is about 1.8 MeV. 

Figure [2] shows the ratios R of the measured (N exp ) and calculated (7V ca i) number of 
electron antineutrino events in the reactor experiments Bugey-4 132|, R0VN091 133 
Bugey-3 134 , Gosgen 135 , ILL 136 , Krasnoyarsk 137 , Rovno88 138 , SRP 139 


Chooz 1401, Palo Verde 141 , Double Chooz [142j, and Daya Bay 122 124 at different 
distances L. As shown in the figure, the average ratio is R — 0.933 ± 0.021, which 
indicates a deficit with a nominal statistical significance of about 3.1cr. 

The reactor antineutrino anomaly can be explained by neutrino oscillations with 
an oscillation length which is shorter than about 20 m. Since the relation between a 
squared-mass difference Am 2 and the corresponding oscillation length L osc is 


L osc = 


AnE 
Am 2 


~ 2.5 


E [MeV] 
Am 2 [eV 2 


m. 


(40) 


and the average energy of the antineutrinos detected in a reactor experiment is about 4 

9 We call “nominal” the statistical significances of the indications in favor of short-baseline neutrino 
oscillations in order to emphasize that they depend on the estimated standard uncertainties. 
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Figure 3. (a) 71 Ga —> 71 Ge transitions induced by 51 Cr and 37 Ar electron neutrinos, 
(b) Ratios R of the measured (iV exp ) and calculated (N ca \) number of electron neutrino 
events in the GALLEX and SAGE radioactive source experiments. The horizontal 
shadowed red band shows the average ratio R and its uncertainty. For each experiment 
the error bar shows the experimental uncertainty. 


MeV, the required squared-mass differences is 

Am| BL £ 0.5 eV 2 . 


(41) 


3.2. The Gallium neutrino anomaly 


106,143-148 


The Gallium neutrino anomaly 
v e measured in the Gallium radioactive source experiments GALLEX |149[jl51 


is a short-baseline disappearance of 

and 


SAGE [l43|p2Hl54 


In these radioactive source experiments the detectors of the GALLEX and SAGE 
solar neutrino experiments have been tested with intense artificial 51 Cr and 37 Ar 
radioactive sources which produce electron neutrinos through the electron captures 


e~ + 51 Cr —► 51 V + i/ e , e“ + 37 Ar 37 C1 + i/ e . (42) 


The neutrino energies and branching ratios are given in Tab. [2} In each experiment 
the radioactive source was placed near the center of the detector and electron neutrinos 
have been detected with the reaction 


v e + (1 Ga —» a Ge + e 


71/ 


(43) 


Figure 3(a) shows the transitions from the ground state of ' 3 Ga to the ground state 
or one of two excited states of (1 Ge which are energetically allowed in the detection 
process. The total detection cross section is given by 


BGT 175 BGTsoo 

O' — CTfrs I 1 + £l75 nnm-^ S500 


BGT 


gs 


BGT 


gs 


( 44 ) 










































Light sterile neutrinos 


14 


Reference 




Method 

bgt 175 

BGT 5 oo 

Krofcheck et a 

T 

157 

158 

71 Ga(p, n) 71 Ge 

< 0.005 

0.011 ±0.002 

Frekers et al. 

15£ 



71 Ga( 3 He, 3 H) 71 Ge 

0.0034 ±0.0026 

0.0176 ±0.0014 


Table 3. Values of the Gamow-Teller strengths of the transitions from the ground 
state of 71 Ga to the two excited states of 71 Ge at 175 keV and 500 keV. 


where <j gs is the cross sections of the transitions from the ground state of a Ga to the 
ground state of a Ge, BGT gs is the corresponding Gamow-Teller strength, and BGTi 75 
and BGT500 are the Gamow-Tellcr strengths of the transitions from the ground state of 


a Ga to the two excited states of ' 4 Ge at about 175 keV and 500 keV (see Fig. 3(a)). 
The coefficients of BGTi 75 /BGT gs and BGT 50 o/BGT gs are determined by phase space: 


6 75 ( 51 Cr) = 0.669, ^5oo( 51 Cr) = 0.220, £i 75 ( 37 Ar) = 0.695, ^oo( 37 Ar) = 0.263 155 . 

The cross sections of the transitions from the ground state of ' 4 Ga to the ground 
state of ,4 Ge have been calculated accurately by Bahcall ||l55 : 


cr g s( 51 Cr) = 55.3 x 10 


- 46 cm 2 , 


cr gs ( 37 Ar) = 66.2 x 10 


- 46 cm 2 , 


and [106,156] 


BGT gs = 0.0871 ± 0.0004. 


(45) 


(46) 


The Gamow-Tcller strengths BGT 175 and BGT 50 o have been measured in 1985 in the 


(p, n ) experiment of Krofcheck et al 
precisely, in 2011 in the ( 3 He, 3 H) experiment of Frekers et al 
listed in Tab. [SI 


157,158] (see Table I of Ref. 155| ) and, more 


159 . The results are 


Figure 3(b) shows the ratios R of the number of electron neutrino events (N exp ) 


measured in the GALLEX and SAGE radioactive source experiments and that calculated 
(A" C ai) with the Frekers et al. Gamow-Teller strengths. The average ratio shown in the 
figure is R = 0.84 ±0.05, which indicates a deficit with a nominal statistical significance 
of about 2.9a. 

Since the average neutrino traveling distances in the GALLEX and SAGE 
radioactive source experiments are (L)gallex = 1-9 m and (L)sage = 0.6 m, from 
Eq. (40) and Tab. [2] one can estimate that the Gallium neutrino anomaly can be 
explained by neutrino oscillations generated by a squared-mass difference 


Am sbl £ 1 eV 


(47) 


3.3. The LSND anomaly 

The LSND experiment 160,161 observed an excess of V e events in a beam of produced 
by /i + decay at rest, 

p + ->• e + + u e ± P M . (48) 
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3+1 

3+1 

3+1 

3+1 

3+2 

3+2 


GLO 

PrGLO 

noMB 

noLSND 

GLO 

PrGLO 

Xmin 

306.0 

276.3 

251.2 

291.3 

299.6 

271.1 

NDF 

268 

262 

230 

264 

264 

258 

GoF 

5% 

26% 

16% 

12% 

7% 

28% 

(Xmin) APP 

98.9 

77.0 

50.9 

91.8 

86.0 

69.6 

(Xmin)DIS 

194.4 

194.4 

194.4 

194.4 

192.9 

192.9 

AxJ-g 

13.0 

5.3 

6.2 

5.3 

20.7 

8.6 

NDFp G 

2 

2 

2 

2 

4 

4 

GoF pg 

0.1% 

7% 

5% 

7% 

0.04% 

7% 

Ay no 

49.2 

47.7 

48.1 

11.4 

55.7 

52.9 

NDFno 

3 

3 

3 

3 

7 

7 

ncr no 

6.4cr 

6.3cr 

6.4cr 

2.6a 

6.1(7 

5.9(7 


Table 4. Results of the fit of short-baseline data taking into account all MiniBooNE 
data (GLO), only the MiniBooNE data above 475 MeV (PrGLO), without MiniBooNE 
data (noMB) and without LSND data (noLSND) in the 3+1 and 3+2 schemes. The 
first three lines give the minimum y 2 (y 2 nin ), the number of degrees of freedom (NDF) 
and the goodness-of-fit (GoF). The following five lines give the quantities relevant for 
the appearance-disappearance (APP-DIS) parameter goodness-of-fit (PG) 164 . The 
last three lines give the difference between the y 2 without short-baseline oscillations 
and Xmin (^Xno)i the corresponding difference of number of degrees of freedom 
(NDFno) and the resulting number of cr’s (nemo) for which the absence of oscillations 
is disfavored. 


The energy spectrum of u fJi is (f>p (E) oc E 2 (3 — 4 E/m^) (see Ref. 
energies E smaller than 


131 


) for neutrino 


E milx = (rrifj, - m e )/2 ~ 52.6 MeV. 


(49) 


The v e events have been detected at a distance L ~ 30 m through the inverse 


neutron decay process (38) in a detector Riled with liquid scintillator in the range 


20 < E e < 60 MeV for the energy E e of the detected positron. 

The nominal statistical significance of the LSND v e appearance signal is of about 


3.8a. However, it must be noted that the similar KARMEN experiment 162,163 did 


not measure any excess of v e events over the background at a distance L ~ 18 m. 


From Eq. (40) one can estimate that the LSND v e appearance signal can be 


explained by z+ —» v e oscillations generated by a squared-mass difference 


+bl £ 0.1 eV . 


(50) 
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Many analyses of short-baseline neutrino oscillation data have been done since the 


discovery of the LSND anomaly in the middle 90’s 80, 81,83, 88, 90, 91,93,165-178 


The activity became more interesting after the discoveries of the Gallium neutrino 


anomaly 106 

145 

148 

179 

183 

104 

106 

107 

184- 

188 . 

The last 


data have been presented in Refs. 109,111 . These analyses take into account the final 
results of the MiniBooNE experiment, which was made in order to check the LSND 
signal with about one order of magnitude larger distance (L) and energy (E), but the 
same order of magnitude for the ratio L/E from which neutrino oscillations depend. 
Unfortunately, the results of the MiniBooNE experiment are ambiguous, because the 
LSND signal was not seen in the neutrino mode —> u e ) (189 and the —$■ u e signal 
observed in 2010 [l90| with the first half of the antineutrino data was not observed in 
the second half of the antineutrino data 191 . Moreover, the MiniBooNE data in both 


neutrino and antineutrino modes show an excess in the low-energy bins which is widely 
considered to be anomalous because it is at odds with neutrino oscillations [103, 104j 


10 


In the following we discuss the results of an updat^] of the global fit of short- 
baseline neutrino oscillation data presented in Ref. El in which the data of the 
following three groups of experiments are considered: 


(A) The % 


—>■ Ve appearance data of the LSND 
E776 [194], KARMEN 


161 


163 , NOMAD 195 , ICARUS 196 


MiniBooNE 
and O 


191 


BNL- 


’ERA 197 


experimental 

(B) The following { v e disappearance data: 1) the data of the Bugey-4 
(l33 , Bugey-3 1134 , Gosgen 135 , ILL [l36|, Krasnoyarsk 137 
SRP 


132 


ROVN091 


Rovno88 138 


139 , Chooz 140 , Palo Verde 141], Double Chooz 142 , and Daya 


Bay 122 124jj reactor antineutrino experiments, with the new theoretical fluxes 
12 97-99] (see Section 3.1); 2) the data of the GALLEX 149-151 and SAGE 


143,152 154] Gallium radioactive source experiments with the statistical method 

considering the recent 71 Ga( 3 He, 3 H) n Ge cross section 


discussed in Ref. 148 


measurement in Ref. 
sin 2 2'd 


159 


(see Section 3.2); 3) the solar neutrino constraint on 

i/ e + 12 C 


12 


15| |106[|20QH202| ; 4) the KARMEN (203||204| and LSND (205 


Ng. s . + e 


scattering data HD, with the method discussed in Ref. 


104 


10 


The interesting possibility of reconciling the low-energy anomalous data with neutrino oscillations 
through energy reconstruction effects proposed in Ref. 


192 


193 


We added the data of several reactor neutrino experiments: 


still needs a detailed study, 
the old results of Rovno88 


138 


and 

SRP 139 , which were omitted in our previous analyses 106 111 but have been considered in the fits 


of other authors 97 109 ; the old results of Chooz 140 

and Daya Bay 


142 


122 


124 


the new results of Double Chooz 

12 For simplicity we do not take into account the correct ana 


and Palo Verde 141 as suggested in Ref. 108 


presented in Ref. 198 
global fits. 


(see also Ref. 199 


ysis of the ICARUS and OPERA data 
because it would not change significantly the results of the 
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sin 2 2fl e|1 


sin 2 2it ee 


sin 2 


Figure 4. Allowed regions in the sin 2 2‘d efJj -Am % 1 , sin 2 2'd ee -Am\ 1 and sin 2 2$^- 
Awij planes obtained in the pragmatic 3+1-PrGLO global fit of short-baseline 
neutrino oscillation data compared with the 3er allowed regions obtained from ^v e 
short-baseline appearance data (APP) and the 3er constraints obtained from*^ short- 
baseline disappearance data ( v e DIS)/i/^ short-baseline disappearance data DIS) 
and the combined short-baseline disappearance data (DIS). The best-fit points of the 
PrGLO and APP fits are indicated by crosses. 


(-) 


(C) The constraints on disappearance obtained from the data of the CDHSW 
experiment [206] , from the analysis [901 of the data of atmospheric neutrino 
oscillation experiment^] from the analysis 103,211 of the MINOS neutral-current 
data |212| and from the analysis of the SciBooNE-MiniBooNE neutrino [213] and 


antineutrino 214 data. 


Table [4] summarizes the statistical results obtained from global fits of the data above 
in the 3+1 and 3+2 schemes (see Section [2] for a global fit in the 3+1+1 scheme see 
Ref. 111] ). In the global (GLO) fits all the MiniBooNE data are considered, including 
the anomalous low-energy bins, which are omitted in the pragmatic PrGLO global 


fits 111 . There is also a 3+1-noMB fit without MiniBooNE data and a 3+1-noLSND 


fit without LSND data. 

From Tab]4[ one can see that in all fits which include the LSND data the absence of 
short-baseline oscillations is nominally disfavored by about 6a, because the improvement 
of the x 2 with short-baseline oscillations is much larger than the number of oscillation 
parameters. 

In all the 3+1 and 3+2 schemes the goodness-of-fit in the GLO analysis is 
significantly worse than that in the PrGLO analysis and the appearance-disappearance 


13 


The analysis of the IceCube data 207-210 , which could give a marginal contribution, have not been 


considered because it is too complicated and subject to large uncertainties. 
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parameter goodness-of-fit is much worse. This result confirms the fact that the 
MiniBooNE low-energy anomaly is incompatible with neutrino oscillations, because it 
would require a small value of and a large value of sin 2 2d e/i 103,104], which are 


excluded by the data of other experiments (see Ref. Ill for further details)]^} Note 
that the appearance-disappearance tension in the 3+2-GLO fit is even worse than that 
in the 3+1-GLO fit, since the Ayp G is so much larger that it cannot be compensated 
by the additional degrees of freedom^ Therefore, we think that it is very likely that 
the MiniBooNE low-energy anomaly has an explanation which is different from neutrino 
oscillations. The cause of the MiniBooNE low-energy excess of ry-like events is going to 


be investigated in the MicroBooNE experiment at Fermilab 216,217 , which is a large 


Liquid Argon Time Projection Chamber (LArTPC) in which electrons and photons can 
be distinguished] (see the review in Ref. 1218j). 


In the following we adopt the “pragmatic approach” advocated in Ref. Ill which 


considers the PrGLO fits, without the anomalous MiniBooNE low-energy bins, as more 
reliable than the GLO fits, which include the anomalous MiniBooNE low-energy bins. 
The 3+2 mixing scheme was considered to be interesting in 2010 when the 


MiniBooNE neutrino 189 and antineutrino 190 data showed a CP-violating tension, 
but this tension almost disappeared in the final MiniBooNE data 191 . In fact, from 


Tabj4] one can see that there is little improvement of the 3+2-PrGLO fit with respect 
to the 3+1-PrGLO fit, in spite of the four additional parameters and the additional 
possibility of CP violation. Moreover, since the p-value obtained by restricting the 3+2 
scheme to 3+1 disfavors the 3+1 scheme only at 1.1a, we think that considering the 
larger complexity of the 3+2 scheme is not justified by the datap] 

Figure |d] shows the allowed regions in the sin 2 2t? eM -Am 2 1 , sin 2 2i? f . P -Am 2 1 and 


Ml 


sin 2 2'd MM -Am| 1 planes obtained in the 3+1-PrGLO fit. These regions are relevant, 
respectively, for { u p —* 'u e appearance, ( u e disappearance and { u\ disappearance searches. 
The corresponding marginal allowed intervals of the oscillation parameters are given in 




(-) 


(-) 


Tab. 5 Figure 4] shows also the region allowed by —> v e appearance data and the 
constraints from ( z+, disappearance and 'v Pi disappearance data. One can see that the 
combined disappearance constraint in the sin 2 2 , d efl -Am\ x plane excludes a large part of 
the region allowed by'rv —> v e appearance data, leading to the well-known appearance- 


disappearance tension 14,101 104,109,188,215 quantified by the parameter goodness- 
of-fit in Tab]4l 

It is interesting to investigate what is the impact of the MiniBooNE experiment on 


14 


One could fit the three anomalous MiniBooNE low-energy bins in a 3+2 scheme 14 by considering 


the appearance data without the ICARUS 196 and OPERA 197 constraints, but the required large 


transition probability is excluded by the disapp eara nce data. 
15 This behavior has been explained in Ref. 215 
Ref. 


It was found also in the analysis presented in 


16 


In 


109 

the 


MiniBooNE mineral-oil Cherenkov detector ^-induced events cannot be distinguished from 


^y-induced events which produce only a visible photon (for example neutral-current 7r° production in 
which only one of the two decay photons is visible). 


1 ' See however the alternative discussion in Ref. 109 . 
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CL 

A m 2 41 

[eV 2 ] 

sin 2 2$ eil 

sin 2 2i3 ee 

sin 2 2'd IHl 

68.27% 

1.57- 

1.72 

0.0011 - 0.0018 

0.085 - 0.13 

0.039 - 0.066 

90.00% 

1.53 - 

1.78 

0.00098 - 0.0020 

0.071-0.15 

0.032 - 0.078 

95.45% 

1.50 - 

1.84 

0.00089 - 0.0021 

0.063 - 0.16 

0.030 - 0.085 

99.00% 

1.24 - 

1.95 

0.00074 - 0.0023 

0.054- 0.18 

0.025 - 0.095 

99.73% 

0.87- 

2.04 

0.00065 - 0.0026 

0.046 - 0.19 

0.021 - 0.12 


Table 5. Marginal allowed intervals of the oscillation parameters obtained in the 
global 3+1-PrGLO fit of short-baseline neutrino oscillation data. 


the global analysis of short-baseline neutrino oscillation data. With this aim, we consider 
two additional 3+1 fits: a 3+1-noMB fit without MiniBooNE data and a 3+1-noLSND 
fit without LSND data. From Tab|4] one can see that the results of the 3+1-noMB 
fit are similar to those of the 3+1-PrGLO fit and the nominal exclusion of the case of 
no-oscillations remains at the level of 6cr. On the other hand, in the 3+1-noLSND 
fit, without LSND data, the nominal exclusion of the case of no-oscillations drops 
dramatically to 2.6cr. In fact, in this case the main indication in favor of short-baseline 
oscillations is given by the reactor and Gallium anomalies which have a similar statistical 
significance 106 . Therefore, it is clear that the LSND experiment is still crucial for the 
indication in favor of short-baseline —y P e transitions and the MiniBooNE experiment 
has been rather inconclusive. 


5. /? decay and neutrinoless double-/? decay 


The existence of massive neutrinos at the eV scale can be probed in /?-decay experiments 


219-225 and in neutrinoless double-/? decay experiments 13,106 


226-234 


5.1. Tritium (3 decay 

The most sensitive experiments on the search of the effects of neutrino masses in /? decay 
use the Tritium decay process 


3 H —>■ 3 He + e + v e . (51) 

Non-zero neutrino masses distort the measurable spectrum of the emitted electron. It 
is convenient to consider the Kurie function (see, for example, Ref. [56]) 

K\T) = {Q- T e ) J2 \Uek\ 2 \j(Q - T e f - m 2 k 0(Q - T e - m k ), (52) 

k 

where T e is the electron kinetic energy, Q = M3 h — M3 n e — m e ~ 18.574 keV is the Q- 
value of the process, and 0 is the Heaviside step function. Considering an experiment 
in which the energy resolution is such that mk -C Q — T e for the three standard light 
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K\T) ~ (Q - T e ) \J {Q — T e ) 2 — <d(Q - T e - m p ) 

+ (Q ~ T e ) J2\ U ^\ 2 \I(Q- T e? * m l ®(Q-T e - m k ), (53) 


k> 4 


with the effective light neutrino mass 


1/2 


™/s = ( Lj I Uet\ 2 ml 


(54) 


K k =1 


The most stringent 95% CL limits on mp have been obtained in the Mainz 235 and 
Troitsk 236 experiments: 


rrip < 2.3eV (95% C.L., Mainz), 
mp < 2.1 eV (95% C.L., Troitsk). 


(55) 

(56) 


The experiment KATRIN 237 , which is under construction and is scheduled to start 


data taking in 2016, will aim to reach a sensitivity of 0.2 eV at 90% C.L. for nip in five 
years of running^ 


The expression (53) of the Kurie function shows that a heavy nonstandard neutrino 


mass m k with k > 4 can be measured by observing a kink of the kinetic energy spectrum 
of the emitted electron at Q — mk below the end point 104,222,241-247 . Recently, the 


Mainz 219 and Troitsk 220,221 collaborations reanalyzed the measured spectra which 


led to the limits (55) and (56) searching for the effect of a heavy neutrino mass m^. The 
negative results led to the upper bounds for the mixing factor | t/ e 4 1 J shown in Fig. [HJ 
These bounds imply that the squared-mass difference that can explain the reactor and 
Gallium anomalies (see Sections 3.1 and 3.2) must be smaller than about 40 eV 2 at 90% 
C.L. 


107 


This is encouraging for future radioactive source experiments which aim at 
a “smoking-gun” measurement of short-baseline oscillations as a function of distance 
and/or energy, because the predicted range for the oscillation length L°^ = AnE / Am^ 
is (1071 

L osc 

6 cm < < 3 m (2cr). (57) 

~ E [MeV] v ’ v ’ 

Some studies were performed to analyze the sensitivity of the KATRIN experiment 


237 to the effects of heavy sterile neutrinos with keV-scale masses 248-250 and light 
eV-scale sterile neutrinos 1222— 225]. According to Ref. |225 , in three years of data 
taking the sensitivity of the KATRIN experiment will cover most of the region of the 
mixing parameters of the sterile neutrino allowed by the reactor and Gallium anomalies 
in the 3+1 framework. 


18 Other future projects are described in the reviews 238 ■ 240 
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Figure 5. (a) Upper 90% C.L. limit on sim(d) = |C/ e 4 1 2 as a function of 

m(v 4 ) 2 = m 4 obtained in the Mainz experiment 219 . (b) Upper 95% C.L. limits 


on [/g 4 = | t/ e 4 1 2 as functions of m 2 obtained with different statistical methods in the 


Troitsk experiment 221 


5.2. Neutrinoless double-/3 decay 

The implications of non-standard mainly sterile massive neutrinos at the eV scale 
for neutrinoless double-/? decay experiments have been studied by several authors 


13.106,176.226-234 


If massive neutrinos are Majorana particles (see the recent reviews in Refs. |77,234j), 
in the case of 3+1 mixing the rate of neutrinoless double-/? decay is proportional to the 
square of the effective Majorana mass 


\mpp | = \\U el \ 2 m 1 + \U e2 \ 2 e ia *m 2 +\U ( 


I 2 +“3 


e31 e d m 3 


+ \U el \ 2 e‘“‘ m ,| 


(58) 


In this expression there are three completely unknown complex phases a? 2 , a 3 , a 4 
which depend on the Majorana phases in the neutrino mixing matrix (see, for example, 
Refs. 


56-58)). These unknown complex phases can generate cancellations between the 


different mass contributions. Figure [ 6 ] shows the range of allowed values of \mpp\ as a 
function of the lightest neutrino mass in the cases of 3 v and 3+1 mixing with Normal 
and Inverted Ordering of the three lightest neutrinos 233 . The 3z/ mixing parameters 


are those in Tab. [I] and the sterile neutrino mixing is that obtained in the pragmatic 
3+1-PrGLO global fit of short-baseline neutrino oscillation data discussed in Section |4j 
From Fig j 6 ] one can see that the presence of an additional massive neutrinos at 
the eV scale can change dramatically the predictions for the possible range of values of 
\mpp\ 13,106,176,226 234 . 
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Lightest mass: m, [eV] 


Lightest mass: m 3 [eV] 


Figure 6. Value of the effective Majorana mass rripp as a function of the lightest 
neutrino mass in the cases of 3v and 3+1 mixing with Normal and Inverted Ordering 
of the three lightest neutrinos 233 . The signs in the legends indicate the signs of 


e la2 , e*“ 3 ; e*“ 4 = ±1 for the cases in which CP is conserved. The intermediate yellow 
regions are allowed only in the case of CP violation. 


6. Cosmology 


In this Section we review the cosmological effects of a massive neutrino at the eV 
scale in the standard ACDM cosmological framework (see also the recent reviews in 
Refs. 16,251 and the book in Ref. [252|). Since such a massive neutrino is mostly 


sterile, as explained in Section [2} its mass is traditionally denoted by the symbol m s in 
cosmological discussions. In this section we use this notation, keeping in mind that its 
real meaning in the 3+1 mixing scheme is m s = m 4 . Moreover, in the discussion of the 
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combined analysis of cosmological data and short-baseline oscillation data we consider 


mi, m 2 , m 3 <C m 4 , so that m s = m 4 ~ \J Am 41 = \/A 


m: 


sbl- 


After introducing in Subsection 6.1 the parameterization of the neutrino density, 


in Subsection (L2 we review the effects of light sterile neutrinos on the observables 
generated in the first part of the evolution of the Universe, when the sterile neutrinos 
were relativistic: the Cosmic Microwave Background (CMB) and the nuclear abundances 


produced by Big Bang Nucleosynthesis (BBN). In Subsection 6.3 we discuss the effects of 


light sterile neutrinos on the formation of Large Scale Structures (LSS), which occurred 


after the sterile neutrinos became non-relativistic. Finally, in Subsection 6 A 
the current cosmological bounds on light sterile neutrinos. 


we review 


6.1. Neutrino parameterization 

It is convenient to parametrize the neutrino contribution to the radiation content in the 
early Universe in terms of an effective number of degrees of freedom N e g, such that the 
total energy density of relativistic species p r is given by 


Pr 


7 / 4 \ 4/3 

1 + ~ ' IT) NeS 


Pi — [1 + 0.22711V e ff] p 7 , 


(59) 


where p 7 is the energy density of photons. One unit of N e g corresponds to the 
contribution of one single family of active neutrinos which were in equilibrium in the 
early Universe and passed through an instantaneous decoupling at a temperature of 
about 1 MeV. The factor 7/8 is for a fermionic degree of freedom, while the factor 
T v /T. 7 = (4/11) 4 / 3 comes from the fact that after neutrino decoupling there is a 
entropy transfer from electrons to photons due to e ± annihilations, which enhances 
the photon temperature, but not the neutrino temperature. Since in the real Universe 
the decoupling of neutrinos was not instantaneous and the neutrinos were not completely 
decoupled when the e ± annihilation occurred, the effective number of active neutrino 
is slightly larger than three: N e g = 3.046 253 . Assuming that all the additional 


contributions to the effective number of relativistic species come from sterile neutrinos, 
their contribution to the total radiation energy density is quantified by AN e g = 
N e g — 3.046, which is given by 


rei 

A N eS = ^ = 


I ^ r4 

815 " 


1 -1 


-t [ dpp 3 f,(p), 


(60) 


where p v is the energy density for one active neutrino species, p v s el is the energy density of 
sterile neutrinos when relativistic, p is the neutrino momentum, f s (p) is the momentum 
distribution and T u = (4/11 ) 1//3 T 7 . 

After they become non-relativistic, neutrinos contribute to the matter energy 
density of the Universe. The sterile neutrino contribution can be parameterized in 
terms of the dimensionless number (see, for example, Ref. |254j) 

UJ S = n s h 2 = —h 2 = — r !h. f dpp 2 fs{p), ( 61 ) 

Pc Pc J 
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where p s is the energy density of non-relativistic sterile neutrinos, p c is the critical 
density and h is the reduced Hubble parameter. Alternatively, oo s can be converted in 


the effective neutrino mass 255 


m'f = 94.1(n s eV. 


(62) 


All the quantities that we introduced depend on the neutrino momentum 
distribution f s {jp). It is important that if the light sterile neutrinos decouple from 
the rest of the plasma when they are still relativistic, f s (p) does not depend on m s . 
It depends only on the production mechanism. The simplest example is one species 
of light sterile neutrinos that are generated by active-sterile oscillations in the early 
Universe [256-261] with the same temperature of active neutrinos. In this case A N e g = 1 
and lu s ~ m s /(94.1eV). 

If the light sterile neutrino thermalizes at a temperature T s = aT U: its momentum 
distribution is given by the standard Fermi-Dirac distribution 


Up) = 


e p/T s _|_ ]_ 


In this case, called “thermal scenario” (TH), from Eqs. (60) and (61) we obtain 


A N eS = a 4 , 


m. 


c o* = a 


94.1 eV’ 


= a 3 m s = A N 3 ^ 4 m s . 


(63) 


(64) 


For a non-thermal sterile neutrino generation, there are several possible 
mechanisms. A popular one is the non-resonant production scenario, also called 
“Dodelson-Widrow scenario” (DW) |262|, which is motivated by early active-sterile 
neutrino oscillations in the limit of zero lepton asymmetry and small mixing angle. 
In this scenario the sterile neutrino has a momentum distribution 


Up) = 


P 


e p / T " + 1 ’ 

where f3 is a normalization factor, leading to 


AlVefj — f3, u s — (3 


m x 


94.1 eV’ 


mf = (3m s = AN eS m s . 


(65) 


( 66 ) 


The DW and the TH models have an exact degeneracy, since they are related by a = /3 1 / 4 
and m™ = m° w ( d 1 / 4 1263,264b 


6.2. Physical effects as radiation in the early universe 

As relativistic components, additional neutrino degrees of freedom change the time of 
matter-radiation equality, whose redshift z eq is given by 


1 T ^eq 


Pm 

Pr 


m 
CO r 


m 1 

Uy 1 + 0.2271A e g ’ 


(67) 
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where p m and p r are the matter and radiation densities, cu, = Ojh 2 , with Qj = pi/p c , 
and we used Eq. (59). 

A shift in the matter-radiation equality affects the position and the shape of 
the acoustic peaks of the CMB (see Ref. (Hi]). At photon decoupling (also called 
“recombination”) the extra radiation component enhances the expansion rate H = a/a, 
where a(t) is the scale factor in the Friedmann-Robertson-Walker metric (see, for 
example, Ref. 252 ). This increase of H generates a decrease of the comoving sound 
horizon r s oc H 265 and a reduction of the angular scale of the acoustic peaks 


6 S = r s /D a, where Da is the angular diameter distance, leading to a shift of the CMB 
peaks towards higher multipoles (see Fig. 2(a) of Ref. [l6|). In addition, if matter- 
radiation equality is delayed, the amplitude of the first CMB peak at t ~ 200 is increased 
by the early Integrated Sachs Wolfe (ISW) effect, since decoupling occurs when matter 
domination is at an earlier stage and the subdominant radiation component causes a 
slow decrease of the gravitational potential (see Figs. 2(a) and 2(b) of Ref. 1161). 

These effects of additional relativistic neutrinos can be partially compensated if 
other cosmological parameters are simultaneously varied. For example, if the total 
matter density co m is also increased without altering the baryon density u;&, so that the 
ratio between odd and even CMB peaks is not altered, according to Eq. (67) z eq can 
be kept fixed and the two effects discussed above do not appear. However, one cannot 
obtain exactly the same CMB spectrum as in the standard case, because additional 


relativistic neutrinos increase the Silk damping effect at high multipoles 265 -267 


Silk damping is a diffusion damping of the oscillations in the plasma which occurs 
at high-multipoles because the decoupling of the baryon-photon interactions is not an 
instantaneous process. During the time when the decoupling occurs, radiation free- 
streams and the temperature fluctuations at scales smaller than the radiation free- 
streaming scale are damped, because on such scales photons can move freely from 
underdensities to overdensities and vice versa. The damping depends on the ratio r^/r s , 
where rq oc iW 1 / 2 is the photon diffusion length at recombination 265 . Since at 
fixed z eq we have H 2 oc p r = (1 + 0.2271A e fj) p 7 , an increase of N e g corresponds to an 


increase of H and an increase of rq/r s oc H 1 ^ 2 , which enhances the Silk damping at 


high-multipoles 265 


Another important effect is related to BBN: the number of relativistic degrees of 
freedom fixes the expansion rate during BBN, that in turn fixes the abundances of 
light elements. BBN can thus give strong constraints on N e g from the observations 

According to Ref. 


of the primordial abundances of light elements 268-271 


272 


BBN limits the effective number of additional relativistic species to AJV e g < 1 at 95% 
C.L., regardless of the inclusion of CMB constraints on the baryon density Vt^h 2 . More 
recently, the authors of Ref. 273 obtained A N e g < 0.2 at 95% C.L. considering the 
BBN and CMB data. 
















Light sterile neutrinos 

6.3. Physical effects as massive component 


26 


Massive neutrinos have an effect in the early Universe due to their velocity. Until they are 
relativistic they free-stream over distances of the order of the Hubble radius if -1 . After 
they become non-relativistic their comoving free-streaming length diminishes together 
with their velocity. As a consequence, there is a maximum comoving free-streaming 
length at the time of the non-relativistic transition, which corresponds today to the 
length A nr = 27r/fc nr , with the wavenumber k nv given by 252 


o.oi78 nil 2 ( y 


1/2 


m s x 1 / 2 
leV7 


h Mpc 


-i 


( 68 ) 


The neutrino free-streaming produces a suppression in the matter power spectrum 
at scales smaller than A nr which has two causes: the absence of the contribution of 
neutrino perturbations and a suppression of the growth of Cold Dark Matter (CDM) 
perturbations (see, for example, Refs. 56,251,2521). 


6.J h Current bounds from cosmology 

The best measurements of the CMB anisotropies come from the second data release of 
the Planck satellite experiment (Planck 2015) |76,274 , which improve the results of the 
first data release (Planck 2013) |255„275j. Concerning neutrinos, in Ref. 1761 the Planck 
collaboration presented constraints on the sum of the active neutrino masses on 

the effective number of relativistic species N eS , plus joint constraints on N e g) 

and on (m® ff , A N ( 


eff 


In the analysis of the Planck collaboration 76 the constraints on N e g alone have 


been obtained considering massless neutrinos in the ACDM +N e g model. However, it is 
interesting to present them in order to understand how the constraints on N e g change 
when different experimental results are taken into account. Considering the temperature 
(Planck TT) and the low-£ polarization (lowP) data, N e g = 3.13 ± 0.32 76 , which is 


consistent with the standard three-neutrino value N e g = 3.046. The inclusion of BAO 


observations 276-278 tightens slightly the constraint to N e g = 3.15 ± 0.23, leading to 


76 


the upper bound A N e g < 1 at more than 3a 

Considering both and AA e ff as free parameters in the fit of cosmological data 


in a ACDM +A r e ff+m(r model, the Planck collaboration obtained 76 


N eS <3.7 


< 0.52 eV 


95%, Planck TT+lowP+lensing+BAO. 


(69) 


In this case, the constraint on N e g is not significantly changed with respect to the ACDM 
-t -N e g model, and A N e g = 1 is still excluded at more than 3cr. High values for A N e g are 
allowed preferably for sterile neutrinos with small effective masses, corresponding to high 
values of the rms amplitude as of linear fluctuations today at a scale of 8 h~ l Mpc. On the 
contrary, for ~ 1 eV one obtains small values of a%, because neutrino free-streaming 
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Figure 7. Allowed la and 2a regions in the mf-N e g plane for the ACDM 
+N e s+m e s B model with three different combinations for the SBL prior on m s : no SBL 
prior (gray), with SBL prior in the DW scenario (red) and in the TH scenario (blue). 
The dotted (dashed) lines correspond to constant values, indicated in eV, for m s in 
the DW (TH) scenario. The green circle represents a fully thermalized sterile neutrino 
with m s = 1 eV. The two figures, adapted from Ref. [291 , have been obtained with 
different cosmological data sets. m CMB data only (Planck 2013 TT at £ < 2500, 
WMAP .E-mode polarization at i < 23, ACT and SPT data at high £). ESI CMB data, 


the prior Hq = 73.8 ± 2.4 Km s 
the 6dFGRS [278 , SDSS-II [293 


-l 


Mpc 1 from HST 
and BOSS DR9 


292 


294 


the BAO measurements of 
surveys, and the cluster mass 
distribution at low and high redshift calculated from the data of the Chandra Cluster 


Cosmology Project 286 with the likelihood implementation of Ref. 287 


suppresses the amplitude of fluctuations at small scales. In fact, as we will discuss in 
the following, a massive sterile neutrino can help solving the current tension between 
CMB measurements and astrophysical determinations of a 8 (the values obtained for <7g 
from weak lcnsing 1279-283 , cluster counts |284,285 , cluster mass functions |286,287 


and redshift space distortions 288-290 are lower than CMB estimates 1761). 

Remember that we are dealing with the effective mass mf r , which is related to m s 


by Eqs. (61) and (62) in the general case, by Eq. (64) in the TH scenario, and by Eq. (66) 
in the DW scenario. Therefore, in spite of the stringent upper bound (69) on mf r , it is 
still possible to have a neutrino with m s ~ 1 eV if AiY e fr is small. However, if the sterile 


neutrinos are generated by active-sterile oscillations in the early Universe 256-261 


they are fully thermalized well before CMB decoupling, resulting in AlV e ff — 1, which 
is disfavored by the bound (69) 19 This problem led several authors to propose new 


19 Previous bounds on A N e g for sterile neutrinos at the eV mass scale have been discussed in 


Refs. 270 295 -299 . 
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mechanisms that can relieve the tension: a large lepton asymmetry 259J, 300-310 
neutrino interactions 


new 


311-321 , entropy production after neutrino decoupling 322 


neutrino decay 323 , very low reheating temperature 324, 325 , time varying dark 


energy components 298 , a larger cosmic expansion rate at the time of sterile neutrino 


production 326 , inflationary freedom 327 


The authors of Refs. [188,215,291,328,329] performed combined fits of cosmological 
data and short-baseline (SBL) oscillation data in the framework of light sterile neutrino 
models. The general result is that the sterile neutrino mass is constrained at the scale of 
1 eV by short-baseline oscillation data and the full thermalization of the sterile neutrinos 


in the early Universe is disfavored by cosmological data. As an example, in Fig. 7(a) 
adapted from Ref. 1293 , we show how the addition of the SBL prior on m s in the DW 
(red) and in the Tff (blue) scenarios changes the allowed regions in the plane 

when the CMB data from Planck 2013 (TT, i < 2500), WMAP (A-mode polarization 
at £ < 23) and ACT/SPT (high-f? temperature spectrum) are considered. The dotted 
(dashed) lines correspond to constant values of m s in the DW (Tff) scenario, with the 
values in eV written near the lines. One can note that since CMB data force the allowed 
regions to be as plotted in gray in Fig. 7(a) N eS and wff cannot be both large at the 
same time. Since the SBL prior requires m s ta leV, in the combined analysis N e g is 
constrained below about 3.75 at 2cr. 

While CMB data disfavor eV-scale neutrino masses and the presence of additional 
light sterile neutrino in cosmology, LSS probes give some hints in favor of massive 
neutrinos. Results in this direction were obtained, for example, in Ref. 284 , suggesting 
a preference for = 0.22 ± 0.09 eV for the active neutrinos if the Planck 2013 

cluster counts through the Sunyaev-Zel’dovich (SZ) effect are considered together with 
CMB and BAO data. Some of the analyses in the 2015 update of Planck SZ cluster 
counts continue to indicate a preference for nonzero neutrino masses (see, for example, 
Fig. 12 in Ref. |285]). These results arise from the fact that the free-streaming of massive 
neutrinos can reduce the tension between the values of <j$ obtained with CMB and LSS 
data 279 -290 . For the light sterile neutrinos, similar analyses were performed in the 
ACDM +N e ff+rri e f model by several authors 291,328-333 , in spite of some tension 


between different datasets 334 


Figure 7(b), adapted from Ref. shows the allowed regions in the rrff-N e s 

plane corresponding to those in Fig. 7(a) with the addition to the cosmological data set 
of the prior on H 0 = 73.8 ± 2.4 Km s^ 1 Mpc -1 from LIST 292 , the BAO measurements 


of the 6dFGRS 278 , SDSS-II 293 , and BOSS DR9 294 surveys, and the cluster mass 
distribution at low and high redshift calculated from the data of the Chandra Cluster 
Cosmology Project 286 with the likelihood implementation of Ref. 287 . Comparing 


the gray allowed regions in Figs. 7(a) and |7(b) one can see that the inclusion of LSS 
data favor values of irnff around 0.2 — 0.6eV. The inclusion of the SBL prior (which 
restricts m s around 1 eV) conserves this indication by restricting the allowed value of 
N e g well below N e g = 4, disfavoring a full thermalization of the sterile neutrino. 

The preference for ~ 0.2 — 0.6 eV is driven mainly by the inclusion of the 
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Project 



neutrino source 

E 

L 

status 





(MeV) 

(m) 


SAGE 

335 

336 

*4 51 Cr 

0.75 

<1 

in preparation 

LENS 

337 


z/ e ,P e 51 Cr, 6 He 

0.75, < 3.5 

<3 

abandoned 


CeLAND 1338,339 


Z4 


144 


Ce 


U 


< 6 


abandoned 


Daya Bay 340 


P e 

144 Ce 

h-i 

bo 

1 

CO 

1.5-8 

proposal 

LENA 

34 i 



V e 

51 Cr, 37 Ar 

0.75, 0.81 

<90 

abandoned 

CeSOX 

342 

-344 

P e 

144 Ce 

1.8-3 

5-12 

in preparation 

CrSOX 

342 

344 

"e 

51 Cr 

0.75 

5-12 

proposal 

JUNO 

345 



P e 

144 Ce 

1.8-3 

<32 

proposal 


Table 6. Main features of proposed source experiments. 


SZ cluster counts from Planck 2013, but similar results are obtained with different 
datasets 291,332. In all these cases, the presence of the light sterile neutrino has 


two effects: its free-streaming reduces the fluctuations at small scales, and hence cr 8 is 
decreased, while the presence of additional radiation in the primordial Universe leads to 
a slight increase of Hq. The results obtained from the LSS measurements, however, can 
suffer from unaccounted systematics or some bias in the determination of the cluster 
masses. Therefore, a definite confirmation or rejection of the LSS preference for a 
massive sterile neutrino must await future observations and analyses. 


7. Conclusions and perspectives 

The reactor, Gallium and LSND anomalies can be explained by neutrino oscillations 
if the standard three-neutrino mixing paradigm is extended with the addition of light 
sterile neutrinos which can give us important information on the new physics beyond 
the Standard Model. 

As discussed in Section [4| the global fits of neutrino oscillation data in the 
framework of mixing schemes with one or more sterile neutrinos suffer from a tension 
between the results of appearance and disappearance short-baseline neutrino oscillation 
experiments. Moreover, as discussed in Section [6j the cosmological data indicate a 
tension between the necessity to have a sterile neutrino mass at the eV scale and the 
expected full thermalization of the sterile neutrinos through active-sterile oscillations 
in the early Universe. Hence, the possible existence of light sterile neutrinos at the eV 
scale is a hot topic of current research and discussions. 

The severe appearance-disappearance tension found in the global analyses of short- 


baseline neutrino oscillation data 109,111 can be alleviated adopting the “pragmatic 
approach” advocated in Ref. Ill , in which the anomalous MiniBooNE low-energy 


excess of u e -like events is neglected. The cause of this excess is going to be investigated 


in the MicroBooNE experiment at Fermilab 216-218 


There is an impressive program of many experimental projects aiming at testing 
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Project 


Nucifer (FRA) J357| 
Stereo (FRA) |12 
Neutrino-4 (RUS) 
Poseidon (RUS) 
DANSS (RUS) 
SoLid (BEL) |361 


358 


359 


360 


362 


362 


Hanbit (KOR) 
Hanaro (KOR) 
Prospect (USA) |363 
CARR (CHN) [364 


Pth 

(MW) 

Af target 

(tons) 

L 

(m) 

Depth 

(rn.w.e.) 

70 

0.8 

7 

13 

57 

1.75 

9-12 

18 

100 

1.5 

6-11 

10 

100 

3 

5-8 

15 

3000 

0.9 

10 - 12 

50 

45-80 

3 

6-8 

10 

2800 

1 

27 

10-23 

30 

0.5 

6 

few 

85 

1,10 

7,18 

few 

60 

~ 1 

7,11 

few 


Table 7. Main features of proposed reactor experiments. 


the short-baseline oscillations due to sterile neutrinos (see also the reviews in Refs. 346- 


351|). It is convenient to divide them in the following three categories (see, however, 


also the proposals in Refs. 352-354 out of these categories): 

Source experiments. These experiments use radioactive sources of v e or F e placed 
near or inside a large detector. Table [6] presents a list of the projects which have 
been proposed (see also Ref. 349] ). To our knowledge, only the SAGE |335, 336 
and CeSOX [342-344] experiments are in preparation. 

Reactor experiments. These experiments use a reactor v e source with a detector 
placed at a distance of the order of 10 m. There are several experiments in 
preparation, as shown by the list in Tab. 0 (see also Ref. [3501). 

Accelerator experiments. These experiments, listed in Tab. [8] (see also Ref. |351|), 
aim at checking the short-baseline —>■ P e LSND signal (see Section 3.3) 

in both neutrino (z/^ —y u e ) and antineutrino (z? M —y u e ) mode (see also the 
nuSTORM proposal in Ref. | 355] ). They will also look for the associated z/ M and Ty 
disappearance (see also the NESSiE proposal in Ref. |356|). 

The aim of most of these experiments is to reveal short-baseline oscillations in a 
robust way by measuring distortions in the neutrino spectrum or variations of the flavor 
neutrino detection probability as a function of distance. In source experiments with 


monochromatic ry’s generated by nuclear electron capture (for example SAGE 335,336 
and CrSOX [342,344|), u e disappearance can be measured as a function of distance 


inside a sufficiently large detector. In source experiments with a continuous D e spectrum 
generated by nuclear (3 decay (for example CeSOX [342-344]) and in reactor experiments 
(for example Stereo 1121) both effects can be measured inside a sufficiently large detector 
with a sufficient energy resolution. Some reactor experiments use two detectors at 
different distances (for example Prospect 363 and CARR |364|) or a movable detector 
(for example DANSS [360]). 
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Project P M target E L 

(MW) (tons) (MeV) (m) 


SBN (USA) [365 










> 0.09 

112,89,476 

~ 800 

110,470,600 

J-PARC MLF 

(JP 

N) 

387 

388 



~ 1 

50 

~ 40 

20 

KPipe (JPN) 

389 









~ 1 

684 

~ 236 

32 - 152 

nuPRISM (JP 

N) 

390 







~ 1 

4000 - 8000 

200 - 1000 

1000 - 2000 

IsoDAR-KamLAlN 

D( 

JPN) 

354 

391 

392 

0.6 

1000 

~ 6.5 

10-40 

IsoDAR-JUNO 

CHN) 

345 

392 


0.6 

20000 

~ 6.5 

20 - 100 

OscSNS (USA) 

393 








1.4 

450 

~ 40 

50-70 


Table 8. Main features of proposed accelerator experiments. 


For accelerator experiments a crucial ingredient for reaching a robust result is the 
presence of “near” and “far” detectors (as, for example, in the SBN |365j experiment). 
The near detector provides a normalization of the neutrino flux and cross section which 
allows to measure the oscillations between the two detectors with small systematic 
uncertainty. 

Since several of the experiments listed in Tabs. [6]-[8] will be performed in the next 
years, there are optimistic hopes that the important open problem of the existence of 
light sterile neutrinos will be solved soon. Let us emphasize that a positive result would 
be a major discovery which would have a profound impact not only on neutrino physics, 
but on our whole view of fundamental physics, because sterile neutrinos are elementary 
particles beyond the Standard Model. Lienee, their existence would prove that there is 
new physics beyond the Standard Model at low-energies and their properties can give 
important information on this new physics. 

For neutrino physics, a discovery of the existence of light sterile neutrinos would 
open a rich field of experimental and theoretical research on the properties of the sterile 
neutrinos, their mixing with the active neutrinos and their role in neutrino experiments 
(e.g. in solar 15,87,106,109 167 200 202,366 367] , long-baseline 84 86 368 369] , and 
atmospheric 83, 166] 208-210,370-375 neutrino experiments) in astrophysics (e.g. in 
supernova neutrino experiments 169,376-383] and indirect dark matter detection |384|), 
high-energy cosmic neutrinos 1257,385,386], and in cosmology (see Section [6]). We think 
that this is a great opportunity which stimulates the great interest in the search for light 
sterile neutrinos. 
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